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Abstract 

Effects of time-delayed feedback on pattern formation are studied both numerically and theo- 
I retically in a bistable reaction-diffusion model. The time-delayed feedback applied to the activator 
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', fronts are changed. If the intensities are different, both the critical point of the NIB bifurcation 



and/or the inhibitor alters the behavior of the Nonequilibrium Ising-Bloch (NIB) bifurcation. If 
the intensities of the feedbacks applied to the two species are identical, only the velocities of Bloch 



and the threshold of stability of front to transverse perturbations are changed. The effect of time- 
delayed feedback on the activator opposes the effect of time-delayed feedback on the inhibitor. 
When the time-delayed feedback is applied individually to one of the species, positive and nega- 
tive feedbacks make the bifurcation point shift to different directions. The time-delayed feedback 
provides a flexible way to control the NIB bifurcation and the pattern formation. 
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I. INTRODUCTION 



Pattern formation has been of great interest in a variety of chemical, physical, and biolog- 
ical contexts.-""- A chemical front (interface) which connects two different states of system, 
such as the excited and recovery states in excitable media, or the two stable states in 
bistable media, plays an essential role on the pattern formation. Many literatures focus on 
the dynamics of the front.-""— Of particular interest is the front controlled by Nonequilib- 
rium Ising-Bloch (NIB) bifurcation in the bistable media. The NIB bifurcation describes 
a pitchfork bifurcation at which a stationary Ising front becomes unstable and a couple of 
counterpropagating Bloch fronts appear. In bistable Ferrocyanide-Iodate-Sulfite reactions, 
spirals, oscillating spots, and labyrinthine patterns have been observed.—"^ The spirals occur 
in the Bloch region beyond the Nonequilibrium Ising-Bloch bifurcation. As a sparse spiral, 
it results from an axisymmetry breaking of a shrinking ring. The oscillating spots appear 
near but before the NIB bifurcation. The labyrinthine pattern originates from transverse 
instability of a chemical front in Ising region. Similar patterns were also observed by Szalai 
and De Kepper.— These patterns observed in bistable Ferrocyanide-Iodate-Sulfite reactions 
can be explained successfully in terms of a NIB bifurcation in a generic FitzHugh-Nagumo 
model.-^^^ 

Controlling the pattern formation is an important issue for the study of self-organization 
phenomena far away from thermodjTiamic equilibrium. Recently, time-delayed feedback, 
firstly presented by Ott et al to control the chaotic behavior of a deterministic system,-- has 
been used to control the pattern formation successfully. It can control the tip trajectories 
of spirals in a light-sensitive Belousov-Zhabotinsky reaction.— A global feedback can either 
stabilize the rigid rotation of a spiral or completely destroy spiral and suppress self-sustained 
activity in a confined domain of excitable medium.— The spontaneous suppression of spiral 
turbulence based on feedback has been studied experimentally in a light-sensitive Belousov- 
Zhabotinsky reaction and numerically in a modified FitzHugh-Nagumo model.— With the 
global feedback one can manipulate the competition between patterns with different sym- 
metries (hexagons and rolls).— In a delayed optical system, resonant Hopf triads lead to 
drifting rhombic and hexagonal patterns.— Near the co dimension-two bifurcation points, 
the time delay can result in a transition between Turing and Hopf instabilities.—'^ 

Most of the studies of the effects of time-delayed feedback on the pattern formation focus 
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on the dynamics of patterns in excitable and oscillatory media.— How about the 
effects of time-delayed feedback on controlling pattern formation in bistable media? In this 
work, we study the role of time-delayed feedback played on controlling the NIB bifurcation 
in a bistable FitzHugh-Nagumo model. We focus on the investigation of the NIB bifurcation 
and the transverse instability of front when applying time delay. The underlying mechanism 
of successful control is analyzed. 

II. BISTABLE MODEL 

The bistable media are described by a FitzHugh-Nagumo model: 

ut = u-u^ -v + V^M + F, (1) 
Vt = e{u — aiv — ao) + SV'^v + G, (2) 

where the time delay is applied with the forms: 

F = gMt-r)-u{t)), (3) 
G = gMt-T)-vit)), (4) 

here, the variables u and v represent the concentrations of the activator and inhibitor, 
respectively, and 6 denotes the ratio of their diffusion coefficients, r indicates the delayed 
time. Qu and are the feedback intensities of variable u and v, respectively. The small 
value e characterizes the time scales of the two variables, where v remains approximately 
constant Vf on the length scale over which u varies. The system described by Eqs.(l) and 
(2) can be either of excitable, Turing-Hopf, or bistable type. In this paper the parameter ai 
is chosen such that the system is bistable. The two stationary and uniform stable states are 
indicated by an up state and a down state {u-,V-), respectively. The parameter ao 

represents the symmetry of the system. In the following, we only consider the case that the 
system is symmetric, i.e. ao=0, (M+,f+)=-(u_,f _). A front connects the two stable states 
smoothly. It can be either traveling (Bloch front) or stationary (Ising) which is determined 
by the control parameters. 

Because pattern formations in bistable media are sensitive to the initial and boundary 
conditions, we adopt fixed initial conditions during the numerical simulations. In the one- 
dimensional case (200 grids, using Euler method), we focus on the travehng wave with an 
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initial condition as shown in Fig.l (a). In the two-dimensional case (200x200 grids, using 
Peaceman-Rachford alternating-direction implicit scheme), we mainly concentrate on the 
spiral wave and labyrinthine pattern with an asymmetrical initial condition indicated in 
Fig. 5 (a). The boundary conditions are taken to be periodic and no-flux in one and two 
dimensions, respectively. The space step is dx=dy=1.0 length unit and the time step is 
dt— 0.05 time unit in both cases. 

III. RESULTS AND DISCUSSION 

A. Front bifurcation in one dimension 

Firstly, we study numerically the front bifurcation in one dimension. Without the time 
delay, the system follows a NIB bifurcation upon decreasing e, which leads to a transition 
from stationary Ising front to a couple of counterpropagating Bloch fronts. Figure 1 (b) 
shows the time evolution of the two Bloch fronts with the parameters deep into Bloch 
region. The traveling wave propagates at constant speed. Now, we apply a time delay with 
the forms in Eqs. (3) and (4), but still keep the control parameters the same as those in 
Fig. 1 (b). When t=0.2, gu—0.9, Qv—O.I, the original Bloch fronts slow down and finally 
stop at some place as shown in Fig. 1 (c). This means that the time delay alters the point 
of NIB bifurcation, which results in a transition from initial Bloch fronts into Ising fronts. 

We obtain numerically the dependence of the front velocity on the parameter e as shown in 
Fig. 2, in order to investigate the front bifurcation. In the absence of time-delayed feedback, 
the NIB bifurcation occurs at £=0.036 and it is a pitchfork bifurcation as indicated by the 
sohd circle in Fig. 2 (a) and (b). Then, we apply positive feedback to the system with the 
delayed time t=0.2 and the identical feedback intensities gu—gv—^-^- From the numerical 
results represented by the empty circle in Fig. 2 (a), it is shown that the opening angle 
of the pitchfork shrinks. The larger the delayed time r is, the smaller the opening angle 
becomes. So, by using time delay one can reduce the front velocity. It is interesting that 
under the feedback with identical intensities the critical point of the NIB bifurcation doesn't 
vary. If we apply the time delay with different feedback intensities the situations become 
much different. On one hand, if the feedback intensity gu is larger than g^, the bifurcation 
point shifts to the left and the opening angle of the pitchfork reduces to some extent as 
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FIG. 1. Time evolution of a front and a back, (a) Initial condition; (b) r=0.0; (c) r=0.2, gu=0-9, 
gv=0.1. The other control parameters are: ai=2.0, 6=2.0, and £=0.03. The time-space span in (c) 
is the same as that in (b). 

shown by the sohd square in Fig. 2 (a). The front speed in this case is smaller than that in 
the absence of time delay. This case can induce a transformation from Bloch fronts to Ising 
fronts as shown in Fig. 1. On the other hand, if the feedback intensity Qu is smaller than g^, 
the bifurcation point shifts to the right, and the opening angle reduce as indicated by the 
empty square in Fig. 2 (a). The empty square line may intersect with the solid circle line 
(without time delay) for small r. On the right of the crossover point the front speed with 
time delay is larger than that without time delay. It is clearly that under the time delay 
with gu<gv, the initial Ising front loses stability to evolve into Bloch front. In general, upon 
increasing the ratio of the feedback intensities gu/Ov, the bifurcation point shifts to the left 
gradually. 

In the negative feedback case, g'„<0, gv<0, the opposite is true. The opening angle 
enlarges, leading to the increasing of the front speed. Upon increasing the ratio of the 
feedback intensities gu/ Qv, the bifurcation point shifts to the right gradually as shown in 
Fig. 2 (b). 

It shows that the effect of time-delayed feedback on the first variable opposes that on 
the second variable. There exists competition between the two feedbacks on controUing the 
NIB bifurcation. If the feedback intensities acting on the two variables are identical, the 
NIB bifurcation point does not affected by the time-delayed feedback as shown above. 
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FIG. 2. Dependence of the front velocity on the parameter e. (a) Positive feedback, sohd circle: 
r=0.0; empty circle: r=0.2, g„=0.9, 5(„=0.9; solid square: t=0.2, g„=0.9, 5(„=0.1; empty square : 
r=0.2, (7u=0.1, 5(^=0.9. (b) Negative feedback, solid circle: r=0.0; empty circle: t=0.2, g(„=— 0.9, 
5t,=— 0.9; solid square: t=0.2, (/t(=— 0.9, (71,=— 0.1; empty square : t=0.2, (7u=— 0.1, 5^=— 0.9. The 
other parameters are: 5=2.0, and ai=2.0. 



If the feedback is applied individually, such that (7^=0, or 5't,=0, we can still realize the 
shift of the critical point of NIB bifurcation. Increasing the feedback ((7^=0), for instance 
from negative to positive values, the bifurcation point shifts from right to left gradually. 
On the contrary, if increasing the feedback (5'u=0) from negative to positive values, the 
bifurcation point shifts from left to right. It shows that the effect of the time delay with 
positive feedback on the variables opposes the effect of time delay with negative feedback 
on the variables. Therefore, by using time delay with appropriate forms one can control the 
front bifurcation efficiently. 

In the absence of the time-delayed feedback the front bifurcation in one dimension is 
determined by the relation between the front velocity and the parameter 
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where q = ^J ai + \- In Eqs. (3) and (4), if the delayed time r is small, we can expand 

u{t — t) and v{t — r) as, 



u{t — t) = u{t) 



v{t — t) — v{t) — T 



— T 



du{t) 

~dr 

dvjt) 
dt 



(7) 



(6) 



So, we can obtain: 



(1 + rgujUf = u — — V + V^w, 
(1 + rgv)vt = e{u — aiv) + SV'^v. 



(8) 
(9) 



It shows that the time delay affects the temporal scales of the variables. We use the singular 
perturbation analysis to study the front bifurcation assuming e/5<^l. Following Eq. (5) it 
is straightforward to obtain the implicit expression of the front velocity: 



Obviously, the feedback parameters r, Qu and affect the behaviors of the front bifurcation. 

Figure 3 plots the dependence of the front velocity on the parameter e based on the Eq. 
(10). In the case of identical feedbacks, such that gu=9v, the velocity of Bloch fronts can 
be rescaled. If gu=gv>0, the front velocity decreases that leading to the reduction of the 
opening angle of the pitchfork as shown in Fig. 3 (a). On the contrary, if gu—gv<0, the final 
front velocity increases, which leads to the increasing of the opening angle as indicated in 



Next, we focus on the critical point of front bifurcation, at which c—0. Thus, Eq. (10) 
can be reduced to: 



c(l + rgu) 



3c(l + rgy) 



(10) 



V2q^[c^{l + Tg^y + AeSq^]^ 



Fig. 3 (b). 
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FIG. 3. Plots of Eq. (10) in the (c, e) plane. The control parameters and the notations correspond 
to those in Fig. 2 

In the absence of the time delay, such that t—0, we denote the critical value of the front 
bifurcation by £co- It can be found from Eq. (11) that if gu>gv, the critical value shifts to the 
left. On the contrary, if gu<gv: it shifts to the right. There exists the competition between 
the feedback gu and g^ on controlling the front bifurcation. The corresponding bifurcation 
diagrams are shown in Fig. 3 (a) and (b), respectively. It can explain well the numerical 
results in Fig. 2. The front bifurcations, both without and with the time delay, are plotted 
in the S — e plane as indicated by the thick dash line Sp and the thick dash dot line 5fd in 
Fig. 7. 

Applying the timc-dclaycd feedback to the variables is equivalent in some sense to chang- 
ing their diffusion coefficients. For example, applying the time-delayed feedback with |1 -|- 
T5f„|<|l -|- T^f^,! is equivalent to slowing the diffusion of the inhibitor [see Eq. (11)], there- 
fore, increasing the velocity of Bloch front. Because the width of Bloch front is in inverse 
proportion to the diffusion coefficient 5, applying the time delay with |1 -|- Tg'u|<|l -|- rg^l 
widens the Bloch front. Figure 4 shows the profiles of Bloch fronts without and with the 
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FIG. 4. Profiles of Blocfi fronts without (a) and witli (b) the time-delayed feedback. The tailing 
of Bloch front Li in (a) is smaller than L2 in (b). The width of up state A2>Ai. The feedback 
parameters in (b) are: t=0.2, gu=0-l, gv=0.9. The other parameters are: 5=2.0, and ai=2.0. 

time-delayed feedback. Obviously, the tailing L2 in (b) is wider than Li in (a) and the Bloch 
front is widen by the time-delayed feedback. Therefore, the width of the up state is widened, 
A2>Ai, resulting in an increase in the wavelength of the spiral wave in two dimension. On 
the contrary, if |1 + r(7u|>|l + rg^l, the opposite is true. 

We want to point out here that with extensive numerical simulation the above results are 
still correct when long delays are applied. The feedbacks gu and g^ can be applied either 
individually or simultaneously, which depends on their values. For example, if individual 
feedback g^ with large delay is applied, r=20, (?„=0.0, and gv=0.1, the bifurcation point 
still shifts right, which is equivalent to the manipulation r=0.2, (7„=0.0, and (7t,=10. The 
larger the product T*g^ is, the farther the bifurcation point deviates from the critical point 
Eco- This is because that the delayed time r and the feedback intensity g^ (gv) are coupled 
together as indicated in the derived Eqs. (8) and (9). Obviously, the action of long delay with 
weak feedback intensity is equivalent to that of short delay with strong feedback intensity. 
This provides guidance in practical application. In order to keep the applicability of Taylor 
expansion on deriving the Eqs. (8) and (9) and the consistency between the analytical and 
numerical results, we used small delays throughout the manuscript. 
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B. Transverse instability of front in two dimensions 

In one dimensional case we concentrate on the front bifurcation by analyzing the relation 
between the velocity of front and the parameter e. A planar front could become curve in two 
dimensions. It is necessary to consider further the stability of a planar front to transverse 
perturbation, i.e. the transverse instability of planar front. In this section, we firstly obtain 
patterns deep into the Bloch and Ising regions, and near the NIB bifurcation point without 
time delay. Then we study the effects of time delay on the transverse instability of front. 
Our emphasis is on controlling the transverse instability of front by applying appropriate 
time delay. 



FIG. 5. Numerical simulations of the symmetrically bistable model, (a) Initial condition; The dash 
line (dot line) represents an interface of variable u (v) , which separates the up state (the left part) 
and the down state (the right part), (b) Spiral wave, £=0.02; (c) Stationary Labyrinthine, £=0.05; 
(d)-(f) are the snapshots of the evolvement of breathing labyrinthine at £=0.036, t=250, 650, 2000 
time units. The other parameters are: 6=2.0, and ai=2.0. Grid size: 200x200 space units. 

In order to illustrate and compare the results clearly, we use the same initial condition as 
shown in figure 5 (a). The intersection point of the contours of u and v servers as an initial 
tip for the spiral formation. From Eq. (10) it can be seen that the parameters e and S are 
coupled together. In the following, for simplicity, we keep the parameter 6 constant. 

Deep into the Bloch region, a couple of Bloch fronts counterpropagate and form spiral 
wave as indicated in Fig. 5 (b). The fronts are stable to transverse perturbations. In the set 
of present parameters the obtained spiral is a dense spiral (the up state and the down state 
are symmetric except an angle separation of tt). Deep into the Ising region, starting from 
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the initial condition [Fig. 5 (a)], the part near the domain center firstly evolves into a spiral 
head. Then, the part behind the spiral head undergoes transverse instability and the fronts 
interplay with each other, which resulting in a stationary labyrinthine pattern finally [Fig. 
5 (c)]. In the stationary labyrinthine pattern the up states keep interconnection and own 
identical widths. This process is similar with the observation in Refs. 15, 20. Near the NIB 
bifurcation, the situation becomes more complex, where we observe a breathing labyrinthine 
pattern. Fig. 5 (d)-(f) show three snapshots of the evolvement of breathing labyrinthine 
at t=250, 650, 2000 time units. In this case, the up states can breakdown and reconnect. 
Together with repulsive interaction between fronts, the widths of the up states increase and 
decrease periodically, leading to the formation of a breathing labyrinthine pattern. The 
most difference between the breathing labyrinthine and the stationary labyrinthine is that 
in breathing labyrinthine case the up state does not interconnect entirely and its width 
changes periodically. The present dynamics is similar with that of oscillatory spots.— 




-1 1 

9u 



FIG. 6. Phase diagram spanned by the feedback intensities gu and gv S-Spiral; BL-Breathing 
Labyrinthine; SL-Stationary Labyrinthine. The other parameters are: 6=2.0, oi=2.0, e=0.036, 
and T=0.2. 

In order to illustrate clearly the transformation between various patterns controlled by 
the time delay, we still use the above individual parameter sets and the initial condition [Fig. 
5 (a)]. Starting from the parameters set in Fig. 5 (f), in which a breathing labyrinthine forms 
in the absence of the time delay, if the feedback intensity gu>gv the given initial condition 
evolves into stationary labyrinthine pattern. However, if gu<gv it transits into the spiral 
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pattern. So, the time delay can alter the critical value of transverse instability of planar 
front (see the analysis below, Fig. 7). When applying the time delay to the system with 
the parameters in Fig. 5 (b) and the same initial condition [Fig. 5 (a)], upon increasing 
the ratio of Qui Qv, it will develop into breathing labyrinthine and stationary labyrinthine 
patterns successively. Similarly, if decreasing the ratio of Qu/ Qv with the parameters as in 
Fig. 5 (c), breathing labyrinthine and spiral patterns form in sequence. Figure 6 shows 
a phase diagram spanned by the feedback intensities Qu and g^, in which the gray region 
represents the breathing labyrinthine pattern. It should be mentioned that the boundary 
between spiral patterns and breathing labyrinthine patterns is not sharp because near this 
boundary the arm of the spiral far away the tip could reflect upon touching the domain 
boundary which leading to the breakdown of the arm. Here, we plot the boundary at which 
perfect spirals could form. The wavelength of spiral can be adjusted by varying the feedback 
parameters, as we have depicted above in the one dimensional case. The spiral period is 
around 160 time units. So, the time delay is still applicable for controlling spiral patterns. 
Therefore, by varying the ratio Qu/ Qv one can realize the control of transverse instability of 
planar front. 

0.4 

0.2 
o 0.0 
-0.2 
-0.4 

-0.1 0.0 0.1 

K 

FIG. 7. Dependence of the front velocity on the curvature. Solid line (0): without time delay. 
Dash dot line (1): r=0.2, 5„=0.1, 5„=0.9; Dotted line (2): r=0.2, £/„=0.9, 5^=0.9; Dash line (3): 
T=0.2, 5u=0.9, g'„=0.1. The other parameters are: 5=2.0, ai=2.0, £=0.036. 

In two dimensions the front velocity is modified by the curvature of front. We should 
consider the transverse instability of planar front beside the NIB bifurcation. Here, we use 
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the algorithm in Ref. 20 to analyze the transverse instability of both Ising and Bloch fronts. 
Under the modification by curvature, Eq. (10) can be written as: 



Cr{l + TQu) + K 



3(cr(l + rg^) + Sk) 



(12) 



1 5 



here, is the normal velocity, and k presents the curvature. Figure 7 shows a velocity- 
curvature relation without and with the time delay. The solid line (0) represents the breath- 
ing labyrinthine pattern at NIB bifurcation without time delay [Fig. 5 (f)]. At the center 
of the plot, the slope of the curve indicates critical stability. If applying time delay with 
identical intensities, such that gu=gv, the velocity changes, but the slope of the curve at the 
center still keeps constant [dotted line (2)]. The stability of front to perturbation hardly 
varies. So, one can still observe breathing labyrinthine pattern. When the feedback in- 
tensities gu>gv, the above slope is positive [dash line (3)]. A front becomes unstable to 
perturbation, and it finally evolves into stationary labyrinthine pattern. On the contrary, if 
9u<gv, the mentioned slope becomes negative [dash dot line (1)], and a front keeps stable 
upon suffering perturbation. We can obtain spiral pattern as shown above. 

We now analyze further the stabilities of both Bloch and Ising fronts to perturbation when 
applying the time delay. If the curvature is small, the normal velocity can be replaced by 
Cr=co—dK, in which cq indicates the velocity of planar front. Here, the reduced parameter 
d is not anymore a simple diffusion coefficient of activator as in excitable system.— Its sign 
determines the stability of a front to transverse perturbations. Inserting into Eq. (12) 
and taking Taylor expansion, we can obtain the implicit expression about d: 



It shows that the reduced parameter d is related with the control parameters 6, e, ai, and 
the feedback parameters in the model. If d is negative, the front becomes unstable upon 
suffering transverse perturbations resulting in the labyrinthine pattern as shown in Fig. 5 
(c). If d is positive, the front keeps stable to transverse perturbations leading to the spiral 
wave as shown in Fig. 5 (b). 

For the Ising front Cq=0, so, we have: 



1 - rf(l + r^„) 



3{6 - d{l + rg,)) 3cl{l + Tg,)^{6 - d{l + rg,)) 



(13) 



V2q^[cl{l + rg.f + Aedq^Y'^ y2g'[cg(l + rg^f + 4e5g2]3/2 " 



l-d{l + rgu) 



3{6 ~d{l + Tg,)) 



(14) 
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FIG. 8. Phase diagram for the NIB bifurcation and the transverse instabiUty. The soUd Une 5i 
indicates the Ising front. The thick dash line dp and dash dot Hne 6fd present the front bifurcations 
without and with time delay, respectively. The thin dash line 5b and dash dot line 8bd present 
the Bloch front without and with time delay, respectively. The feedback intensity in (a) 5^=0.9, 
5^=0.1; (b) 5u=0.1, g'„=0.9. The other parameters are: ai=2.0, and t=0.2. 



At the critical point d=Q, the Ising front will undergo transverse instability. So we obtain 
the critical line for the Ising front: 

Si = (15) 

It can be seen that the transverse instability boundary for the Ising front is unaffected 
by the time delay as shown by the solid hne Sj in Fig. 8. For the Bloch front, cqt^O. At 
the critical point to transverse perturbation (d—O), we obtain an implicit expression for the 
Bloch front: 



^(I^Sbdt-Z + T~Z T~Z ^ 1- 
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The positive solution of Sbd defines a boundary of tlie transverse instability of Bloch 
front as shown by the thin dash dot lines in Fig. 8. From Eq. (16) it is found that the 
competition between and alters the boundary. If gu>gv {.9u<9v), the boundary moves 
down (up) as shown in Fig. 8 (a) [Fig. 8 (b)]. When gu=gv the boundary stays constant as 
the case without the time delay, which means that the time delay do not affect the critical 
stability of Bloch front to transverse perturbations if the feedback intensity gu equals to g^. 

IV. CONCLUSION AND REMARKS 

In this work, we have studied the effects of the time-delayed feedback on the NIB bifurca- 
tion in a bistable medium. The results have shown that the time-delayed feedback applied 
to the activator and/or the inhibitor changes the critical point of NIB bifurcation. The 
time delay alters the temporal scales of the reactions, therefore the velocity of Bloch front. 
Large delay with weak feedback intensity is equivalent to small delay with strong feedback 
intensity. The effect of time-delayed feedback on the activator opposes that on the inhibitor. 
So there exists competition between the two feedbacks on controlling the NIB bifurcation. 
Upon increasing the ratio gu/ dv, the critical point of NIB bifurcation shifts left which could 
result in a transition from Bloch front to Ising front, and vice versa. When time-delayed 
feedback is applied individually to one of the species, positive and negative feedback make 
the bifurcation point shift to different directions. In the two-dimensional case, the time de- 
lay can change the stability of front to transverse perturbations. If gu<gv, it could stabilize 
the front upon suffering transverse perturbation, and vice versa. In some sense applying the 
time-delayed feedback to species is equivalent to changing their diffusion coefficients. Thus, 
the wavelength of patterns can be controlled by properly using feedback parameters. 

Although this FitzHugh-Nagumo model is a generic model, it has described successfully 
the dynamics of pattern formation in bistable Ferrocyanide-Iodate-Sulfite reactions, such as 
the bistable spirals, oscillating spots, and labyrinthine patternsi^"— . These phenomena have 
been attributed to the NIB front bifurcation. In this paper, we focus on the generalized 
controlling scheme to the NIB bifurcation by applying time-delayed feedback to one or two 
of the variables. The results have shown the flexibility of this strategy on controlling the 
NIB bifurcation, therefore the transformation of patterns. 

Many real chemical experiments, such as the ferroin-, Ru(bpy)3-, and cerium-catalyzed 
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Belousov-Zhabotinsky systems, are sensitive to visible and/or ultraviolet light.— ii^>22iiMiM-i2E 
People have realized controlling of pattern formation by the time-delayed feedback in light- 
sensitive chemical reactions. For example, by projecting the delayed image uniformly from 
the feedback loop to the gel in the Petri dish, Kheowan and Zykov realized the controlling 
of spiral waves in a thin layer of the light-sensitive Belousov-Zhabotinsky reaction.— The 
radius of the attractor for meandering spiral waves can be effectively manipulated by varying 
the delayed time in the feedback loop. Karl Vanag et al observed oscillatory cluster pat- 
terns in a light-sensitive Ru(bpy)3-catalyzed Belousov-Zhabotinsky reaction.— The catalyst 
Ru(bpy)3 is light-sensitive. Thus, a proper illumination of the active chemical substrate 
can be used for spatial control of the inhibiting process (Br~). Our results have also con- 
firmed that applying the time-delayed feedback only to the inhibitor is enough to control 
the pattern formation. In a light-sensitive ferrocyanide-iodate-sulphite reaction, Lee et al 
observed the pattern transformation via NIB bifurcation by changing the flow rate or the 
input ferrocyanide concentration.— Our results have shown that applying the time-delayed 
feedback for controlling the NIB bifurcation, from the experimental viewpoint, is equivalent 
to changing the residence time or the input ferrocyanide concentration. We hope that our 
results can be verified in one of the light-sensitive reactions with patterned (not uniform) 
illumination after feedback loop. The feedback loop should mainly include: 1) CCD camera, 
2) video recorder, 3) computer which implements the algorithm of eqs. (3) and (4) and out- 
puts the results (patterned images with appropriate intensity) to a projector, 4) projector 
which projects the patterned images inputted from the computer to the chemical substrate. 
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